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Abstract 



We extend the Golec-Biernat-Wiisthoff model for virtual photon - 
proton scattering to include the resolved photon component explicitly. 
The parameters of the resolved photon component are taken from the 
, literature, while the parameters of the dipole-nucleon interaction are 



fitted to the HERA data in a selected limited range of x and Q . A 



O ' good agreement with experimental data is obtained far beyond the 

. region of the fit. 

O 

cn 
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1 Introduction 

p • The recent decade of investigating deep inelastic scattering at very small 

1^ , Bjorken x at HERA has provided precise data for the F2 structure function 

or equivalently for cr/o/ at large center-of-mass energies. Many phenomeno- 
logical analyses have been performed in order to fit the data. The theoretical 
^ ' analyses can be divided into two general classes. One group of models tries 

c3 ■ to fit the data using the so-called dipole representation. In this approach, 

initiated by Nikolaev and Zakharov P], one fits parameters of the dipole- 
nucleon interaction [21 El E] as a function of the transverse quark- ant iquark 
distance. Another group of models uses rather the momentum representa- 
tion OinilZI- Still another approach jSj tries to fit the so-called unintegrated 
gluon distributions to the HERA data (see also [Til]). 

The fits in the dipole representation take into account only a simple quark- 
antiquark Fock component of the photon. However, the higher Fock com- 
ponents seem to be important to understand the diffraction [12 in detail. 
The importance of the higher Fock states is at present not fully understood. 
The first theoretical step in going beyond the qq component has been un- 
dertaken only recently However, no quantitative estimates exist up to 
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now. Only a very schematic QCD-inspired model was considered in Ref . |14j . 
On the phenomenological side, the jet production in virtual-photon-proton 
scattering, especially at small photon virtuality, shows clearly the presence 
of the resolved photon component (see e.g. fTH]) which seems impossible 
to be explained with the quark-antiquark component only. The ratio of 
the dijet cross section with x^^^ < 0.75 (resolved component) to that with 
j,OBS ^ Q 'jTg (direct component) has been found to increase as decreases. 
The variable x^^^ is to be interpreted as the fractional momentum of the 
photon taking part in the dijet production. At large photon virtuality the 
resolved photon component disappears. The observed dependence of the 
resolved photon component is roughly consistent with the naive VDM form 
factor llj. The present NLO calculations of jet and dijet production in- 
clude these phenomenological form factors when going from real to virtual 
photons (see e.g. [IE!)- Such a phenomenological factor is then prescribed 
to the structure of virtual photon, and more precisely to the parton distri- 
butions in the virtual photon. The resolved photon component seems also 
crucial for understanding the world data for the F^a;, Q^) — F2{x,Q'^) p. 
All these arguments put into question the simple fits to the total photon- 
nucleon cross section with the colour dipole component alone, and call for a 
multi-component parametrization. 

In the light of the extremely successful phenomenological description in 
Ref. 12] it is interesting to see if any phenomenological two-component model 
can do a better job. It is the aim of this note to analyze phenomenologically 
if such a two-component model can satisfactorily describe the HERA data 
for the photon-proton total cross section. In our exploratory analysis, the 
higher Fock components are parametrized by the standard vector dominance 
cross section, i.e. our model is similar in spirit to the Badelek-Kwiecihski 
model [m. 

2 Formulation of the model 

It is known that the LO total 7*A^ cross section in the so-called dipole or 
mixed representation can be written in the form 



In this paper we take the so-called quark-antiquark photon wave function of 
perturbative form pp. As usual, in order to correct the photon wave func- 
tion for large dipole sizes (nonperturbative region) we introduce an effective 
quark/antiquark mass (mg// = mo). 
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The dipole representation has been used in recent years to fit the 
virtual photon - nucleon total cross section [21 El- The best fit has been 
achieved in the saturation model of Golec-Biernat-Wiisthoff j2]. In their 
approach the dipole-nucleon cross section was parametrized as 



a{x,p) = (To 



1 — exp 



(2) 



where the Bjorken x dependent radius Rq is given by 



\ V2 

^o(x) = -^(- . (3) 



IGeV \Xo J 

Model parameters (normalization constant ctq and parameters xq and A) have 
been determined by the fit to the inclusive data on F2 for x < 0.01 

In the GBW approach, the dipole-nucleon cross section is parametrized 
as a function of Bjorken x. As discussed in it would be useful to have 
rather a parametrization in the gluon longitudinal momentum fraction Xg 7^ 
X instead of the Bjorken x. Then one could use the unintegrated gluon 
distribution which is related to the dipole-nucleon cross section as 

An f d'^Kt / ^ -^^ 2\ 

a(gg)N{Xg,p) = — / [1 - {% p) \ [X g , K^) 

•-> J f^t 

^ ^'^^ [1 - Jo(M] asJ^ixg, kI) . (4) 
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Therefore, we find it more appropriate to parametrize the dipole-nucleon 
cross section as Xg instead of the Bjorken x. This involves the following 
replacement in Eq.Q 

(j(x,p) -> cr(Xg,p) , (5) 

which means also a replacement of x by Xg in Eq. (jH]) • As discussed in ^HI , an 
exact calculation of Xg in the dipole representation is, however, not possible, 
and we approximate Xg {M'^^ + Q'^) / {W'^ ^-Q"^), where M^^ = my{z(l — z)) 
with niq being effective quark mass = mo for u/u and d/d (anti)quarks 
and iTLq = mo+ 0.15 GeV for s/s (anti)quarks. This means that in our 
approach the dipole-nucleon cross section <J{qq)N = o'iW, Q^, z, p). This must 
be contrasted to the approach of Ref. [2], where there is no z dependence of 

Including higher Fock components of the (virtual) photon, one could write 
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somewhat schematically: 

q 

+ 



(6) 



The differentials above dVL2 and dVl^ represent phase space volumes for the 
qq and qqg components, respectively, and to'2, OJ^, represent the corresponding 
sets of kinematical variables necessary to describe the relevant process. The 
second and all subsequent terms are of the type of a resolved photon. A 
rigorous approach to the problem is rather difficult ^Hj and has not been 
pursued numerically. 

We shall not try to follow the theoretical path sketched above. Our aim 
here is somewhat different. We intend to construct a simple two-component 
model. One component of our phenomenological model is the dipole qq com- 
ponent, while the other one is meant to represent the resolved photon ex- 
plicitly. Trying to keep our model as simple as possible and inspired by the 
phenomenological results mentioned in the introduction, we wish to check 
if the standard vector dominance model (VDM) contribution can be a rea- 
sonable representation of the resolved photon. Our approach should not be 
understood as a replacement of the missing terms in Eq.®. In our opinion, 
the VDM contribution under consideration contains nonperturbative terms 
which cannot be easily generated by the formal expansion (jHl). However, as 
already mentioned in the introduction, in many exclusive processes the VDM 
contribution represents phenomenologically the resolved photon fairly well. 

The cross section for the VDM component, equivalently called here the 
resolved photon component, is calculated in the standard way 



a. 



We take the simplest diagonal version of VDM with p, uj and mesons 
included. As discussed recently in the contributions of higher vector 
states are expected to be damped. Above the meson-nucleon resonances it 
is reasonable to approximate 

1 



pN _ cuN _ _ 
'tot — ^tot ~ 2 



TT^P I TV p 

(^tot + ^tot 



with a similar expression for cr^"* j221- A simple Regge parametrization of 
the experimental pion-nucleon cross section by Donnachie and Landshoff is 
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used ^21j. As in [22], we take 7's calculated from the leptonic decays of 
vector mesons, including finite- width corrections. The factor (1-x) is meant 
to extend the VDM contribution towards larger Bjorken x. 



3 Fit to the HERA data 

In the previous section we presented formulae for the virtual photon - nucleon 
cross section. The relation between a]^l^ and F2 is a matter of convention. 
In the so-called Hand convention one obtains 

(^'<5) = ^,^^^l + ^^J-F2(a:,Q). (9) 

If the Gilman convention is used instead, then 



-7o7(w^,Q^) = ^^yi + ■ P^M') ■ (10) 



We transform the structure function data from in the standard but ap- 
proximate way ^ 

arf{W,Q') = ^^^-F,{x,Q^). (11) 

Then we perform two independent fits to the HERA data. In fit 1, only 
dipole nucleon interaction is included (see Eq.((H)) 

FITl: a7,7 = a];^. (12) 

In fit 2 in addition we include the resolved photon component in the spirit 
of the vector meson dominance model (see Eq.(|7))) 

FIT2: a]:r = aJl"" + 4dm ■ (13) 

In these fits we limit to 0.15 GeV^ < < 10 GeV^. The upper limit 
is dictated by the simplicity of our model. It is known that at large photon 
virtualities one has to include QCD evolution |24j, which is ignored in the 
present analysis for simplicity. The maximal Bjorken x in the data sample 
included in our fit is 0.021, and minimal W = 17.4 GeV. With the selection 
criterion specified, we select 159 experimental data points. 



^Both prescription © and converge to the standard formula below in the limit of 
small X. 
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Table 1: Compilation of fit parameters. 



fit 


mo (GeV) 


ctq (mb) 


Xo 


A 


x' 


FITl 


0.10 


17.0 


9.50e-4 


0.302 


8.125 


dipole only 


0.15 


23.5 


2.00e-4 


0.268 


4.764 




0.20 


36.0 


1.95e-5 


0.235 


3.080 


FIT2 


0.10 


7.5 


0.0238 


0.3160 


1.696 


dipole + VDM 


0.15 


8.0 


0.0194 


0.3107 


1.553 




0.20 


8.0 


0.0198 


0.3213 


1.812 




0.30 


15.0 


1.67e-3 


0.250 


1.412 




0.40 


24.0 


2.20e-4 


0.230 


1.505 




0.60 


55.0 


l.lOe-5 


0.230 


4.632 



In Table 1 we present the model parameters obtained from the fit. The 
region of small is sensitive to the value of the effective quark mass. This 
nonperturbative parameter is related e.g. to the confinement and cannot be 
obtained from the first principles. Therefore, in the table we show results 
with different values of this parameter. From the fit we find cj^**^ ^ (^o**^ 
for the same value of the effective quark mass. We have extended the range 
of effective quark masses in fit 2 (dipole+VDM). A good quality fit can be 
obtained in the broad range of mo. The criterion by itself does not allow 
one to answer the question which set of parameters is better. The value of 
the psr degree of freedom is shown in the last column. The value of in 
fit 2 (dipole+VDM) is much smaller than that in fit 1 (dipole only), which 
is not acceptable statistically. This can be taken as the evidence for resolved 
photon component. 

In order to illustrate how well the model parameters can be determined 
from the fit to experimental data, in FigQ we show two-dimensional maps 
of in both cases. Here mo = 0.15 GeV and 0.20 GeV was taken for fitl 
and fit2, respectively. Well defined minima are clearly seen. It can be seen 
from Table 1 and Fig^that the parameter xo changes dramatically when the 
VDM component is included, while A stays almost the same. 

The quality of the fit can be judged by inspecting Figs|21-El Since there 
is a rather weak dependence of the cross section on W, therefore in the fig- 
ures showing dependence both theoretical curves and experimental points 
are rescaled by an extra factor 2"', where n counts the subsequent subsets 
of data of a given W shown in FigOl and El Only the cross sections for the 
lowest energy chosen (W = 18 GeV) are left unchanged. By careful inspec- 
tion of the figures one can see the superiority of the fit 2. In presenting the 
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results we have made an arbitrary choice of niQ. The results for other sets 
of parameters (different itlq) are almost indistinguishable in the range of the 
fit. They differ somewhat, however, outside the range of the fit where no 
experimental data are available. The theoretical curves with dipole compo- 
nent only underestimate somewhat the low data. We wish to stress that 
the quality of our fit 1 is worse than that of the original saturation model of 
Golec-Biernat-Wiisthoff [21 . We conclude, therefore, that parameterizing the 
dipole-nucleon cross section as a function of the Bjorken x, instead of Xg, is 
essential for the good quality of the fit in Ref.|2]. 

Having shown that a good-quality two-component fit to the HERA data 
with very small number of parameters is possible, we wish to show a decom- 
position of the cross section into the two model components. In FiglHl and 
Figl7|we show separate contributions of both components as a function of W- 
and Q^, respectively. While at low energy the VDM contribution dominates 
due to the subleading reggeon exchange, at higher energies they are of com- 
parable size. The VDM contribution, being a higher twist effect, dominates 
at small values of photon virtualities. At larger the dipole component 
becomes dominant. This effect is almost independent of energy. 

Up to now we have concentrated on very low-x region relevant for DIS 
at HERA. It is interesting to check what happens if we go to somewhat 
larger Bjorken x or smaller energies W. In this region one cannot neglect the 
valence quark contribution to the cross section. Then the cross section is a 
sum of three components: 

^L^iW, Q') = aj:;{W, Q') + a'yZ,{W, Q') + Q') , (14) 

where the last component is calculated according to Eq. ljllj) with 

F^ix, ^ F2(x, = ^^^^ (^^xuU^, Q^) + ^xd,ai{x, Q')^ . (15) 

We freeze below = 0.25 GeV. In the present calculation we take the 
leading order valence quark distributions from Ref.|2S|. The Q^-dependent 
factor in front of the r.h.s. of Eq. ()15|) is necessary when extrapolating the 
quark contribution to the non-DIS, low-Q^ region (see e.g. [221 )• The param- 
eter Qq (=0.8 GeV^) is taken from a global analysis of the experimental data 
in ^22j. In Fig. IHlwe compare predictions of our two models (fits) also with 
fixed target data j2Zll2H|- The fixed target data are represented by solid sym- 
bols, while the HERA data by open circles. Formula (fTTj) is used to calculate 
both experimental and theoretical cx^jj cross sections. A better agreement is 
obtained with model 2 (dipole-|-VDM), especially at ~ 3 - 5 GeV^, i.e. 
for the NMC data. An overestimation of model 2 at small energies and small 
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photon virtualities may be caused by neglecting a form factor responsible 
for correcting the VDM contribution for finite times of hadronic fluctuations 
|22j . Summarizing, there is a phenomenological evidence for the presence of 
the resolved photon component from the analysis of experimental data for 
F2 in consistency with exclusive reactions. 

4 Conclusions 

Recent fits to the total 7*p cross section in the literature include only the 
quark-antiquark component in the Fock decomposition of the photon wave 
function. The contribution of the higher Fock components, neglected so far, 
is not known and difficult to calculate consistently within quantum chro- 
modynamics. The first trials to include the qqg Fock component within 
perturbative QCD have not be quantified in the literature. Nonperturbative 
effects, not easy to implement within the framework mentioned, can be also 
expected. It is known from the phenomenology of the exclusive reactions that 
the traditional vector dominance model in many cases gives a good estimate 
of the effects characteristic for resolved photon. In this note we have ana- 
lyzed if a two-component model, which includes the qq component and the 
more complicated components replaced by the standard VDM, can provide 
a good description of the HERA data for 7*p scattering. 

In order to quantify the effect of the resolved photon we have performed 
two different fits to the HERA data. In fit No.l we include only the dipole 
component. Here we have used the fiexible and successful parametrization 
of Golec-Biernat and Wiisthoff. In comparison to their fit, in our fit we 
parametrize the dipole-nucleon cross section in terms of a variable which is 
closer to the gluon longitudinal momentum fraction Xg than to the Bjorken 
X. Such a fit is useful on its own, as the corresponding unintegrated gluon 
distribution can be used to estimate cross sections for many exclusive pro- 
cesses. In fit No. 2, in addition we include the VDM component while keeping 
the same functional form of parametrization for the dipole-nucleon interac- 
tion. A better fit is obtained if the resolved photon component of the type of 
VDM is included. When going to slightly larger Bjorken x, the model must 
be supplemented for valence quark contribution. If this is done, the model 
describes also the fixed target data quite well. The two models give different 
predictions in the regions of the phase space where no experimental data are 
available. 

Our phenomenological analysis is only a first step towards a better under- 
standing of the role which the higher Fock components of the photon play in 
both inclusive and exclusive processes. The relation of the phenomenological 
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VDM contribution to the formal expansion discussed in the paper requires 
further study. A numerical calculation of higher-order pQCD effects is called 
for to start addressing this question quantitatively. 
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Figure 1: Two-dimensional maps of P^f degree of freedom for fit 1 (left 
panel) and fit 2 (right panel). Please note a different range of Xq for fit 1 
and fit 2. 
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Figure 2: Quality of fit 1 (qq dipole only) - cross sections as a function of 
W. Lines and sets of experimental data are marked by the value of photon 
virtuality in GeV^. The HERA data taken from f^ . 
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ure 3: Quality of fit 1 (qq dipole only) 
The HERA data taken from W^ . 



cross sections as a function of 
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Figure 4: Quality of fit 2 (qq dipole and VDM) - cross sections as a function 
of W. Lines and sets of experimental data are marked by the value of photon 
virtuality in GeV^. The HERA data taken from 
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Figure 5: Quality of fit 2 (qq dipole and VDM) - cross sections as a function 
of Q^. The HERA data taken from ^26]. 
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Figure 6: Decomposition of total 7*p cross section into dipole (dotted) and 
VDM (dashed) contributions for 4 different values of photon virtuality in 
GeV. 
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Figure 7: Decomposition of total 7*p cross section into dipole (dotted) and 
VDM (dashed) contributions for 2 different energies W. 
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Figure 8: Extrapolation of the models towards fixed target data including 
valence quarks. Solid line represent a sum of the model 1 (left panel) or 
model 2 (right panel) and valence quark contribution. The dashed lines show 
model 1 and model 2 separately. The NMC data are shown by solid circles, 
while the E665 data by solid triangles. The HERA data (open circles) are 
shown for reference. 
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